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ABSTRACT 

We place tight constraints on the growth index 7 by using the recent growth 
history results of 2dFGRS, SDSS-LRG, VIMOS-VLT deep Survey (VVDS) 
and WiggleZ datasets. In particular, we investigate several parametrizations 
of the growth index j(z), by comparing their cosmological evolution using 
observational growth rate data at different rcdshifts. Utilizing a standard like- 
lihood analysis we find that the use of the combined growth data provided 
by the 2dFGRS, SDSS-LRG, VVDS and WiggleZ galaxy surveys, puts the 
most stringent constraints on the value of the growth index. As an example, 
assuming a constant growth index we obtain that 7 = 0.602±0.055 for the con- 
cordance ACDM expansion model. Concerning the Dvali-Gabadadze-Porrati 
gravity model, we find 7 = 0.503 ± 0.06 which is lower, and almost 3er away, 
from the theoretically predicted value of jdgp — 11/16. Finally, based on a 
time varying growth index we also confirm that the combined growth data 
disfavor the DGP gravity. 

Keywords: cosmology: cosmological parameters 



1 INTRODUCTION 



Recent studies in observational cosmology, using all the 
available high quality cosmological data (Type la super- 
novae, cosmic microwave background, baryonic acous- 
tic oscillations, etc), converge to an emerging "standard 
model". This cosmological model is spatially flat with 
a cosmic dark sector usually formed by cold dark mat- 
ter and some sort of dark energy, associated with large 
negative pressure, in order to explain the observed ac- 
celerating expansion of the Universe (cf. Tegmark et al. 
2004; Spergel et al. 2007; Davis et al. 2007; Kowalski et 
al. 2008; Hicken et al. 2009; Komatsu et al. 2009; Hin- 
shaw et al. 2009; Lima & Alcaniz 2000; Jesus & Cunha 
2009; Basilakos & Plionis 2010; Komatsu et al. 2011 and 
references therein). Despite the mounting observational 
evidence on the existence of the dark energy compo- 
nent in the universe, its nature and fundamental ori- 
gin remains an intriguing enigma challenging the very 
foundations of theoretical physics. Indeed, during the 
last decade there has been an intense theoretical debate 
among cosmologists regarding the nature of the exotic 
"dark energy" . The absence of a fundamental physical 
theory, concerning the mechanism inducing the cosmic 
acceleration, has opened a window to a plethora of al- 
ternative cosmological scenarios. Most are based either 
on the existence of new fields in nature (dark energy) 
or in some modification of Einstein's general relativity, 
with the present accelerating stage appearing as a sort of 
geometric effect (for reviews see Copeland, Sami & Tsu- 



jikawa 2006; Caldwell & Kamionkowski 2009; Amendola 
& Tsujikawa 2010 and references therein). 

In order to test the validity of general relativity 
on cosmological scales, it has been proposed that mea- 
suring the so called growth index, 7, could provide an 
efficient way to discriminate between scalar field dark 
energy (hereafter DE) models which admit to general 
relativity and modified gravity models (cf. Ferreira & 
Skordis 2010 and references therein). Linder & Cahn 
(2007) have shown that there is only a weak depen- 
dence of 7 on the equation of state parameter w(z), im- 
plying that one can separate the background expansion 
history, H(z), constrained by a large body of cosmo- 
logical data (SNIa, BAO, CMB etc), from the fluctu- 
ation growth history, given by 7. In this framework, 
it was theoretically found that for those DE models 
which adhere to general relativity the growth index 7 
is well approximated by 7 ~ 1™ Z 5 (see Silveira & 
Waga 1994; Wang & Steinhardt 1998; Linder 2004; Lin- 
der & Cahn 2007; Nesseris & Perivolaropoulos 2008; 
Lee & Kin-Wang 2010), which reduces to 7A ~ 6/11 
for the traditional ACDM cosmology w(z) = — 1. On 
the other hand, in the case of the braneworld model of 
Dvali, Gabadadze and Porrati (2000; hereafter DGP) 
the growth index becomes 7dgp — 11/16 (see also Lin- 
der 2004; Linder et al 2007; Gong 2008; Wei 2008; Fu, 
Wu & Hu 2009), while for the f(R) gravity models we 
have 7 ~ 0.41 — 0.212 for Q, m o = 0.27 (Gannouji, Moraes 
& Polarski 2009; Tsujikawa et al. 2009; Motohashi, 
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Starobinsky & Yokoyama 2010). Indirect methods to 
determine 7 have also been proposed (mostly using a 
constant 7), based either on the observed growth rate of 
clustering (Nesseris & Perivolaropoulos 2008; Guzzo et 
al. 2008; Di Porto & Amendola 2008; Gong 2008; Dos- 
set et al. 2010; Samushia, Percival & Racanelli 2012; 
Hudson & Turnbull 2012) providing a wide range of 7 
values 7 = (0.58 — 0.67)i 011 1 ' 17 , or on massive galaxy 
clusters Vikhlinin et al. (2009) and Rapetti et al. (2010) 
with the latter study providing 7 = 0.42+Q'^g, or even 
on the weak gravitational lensing Daniel et al. (2010). 
Gaztanaga et al. (2012) performed a cross-correlation 
analysis between probes of weak gravitational lensing 
and redshift space distortions and found no evidence 
for deviations from general relativity. With the next 
generation of surveys, based on Euclid and BigBOSS, 
we will be able to put strong constraints on 7 (see for 
example Linder 2011; Belloso, Garcia-Bellido & Sapone 
2011; Di Porto, Amendola & Branchini 2012 and refer- 
ences therein) and thus to test the validity of general 
relativity on extragalactic scales. 

The scope of the present study is along the same 
lines, ie., to place constraints on the growth index us- 
ing a single cosmologically relevant experiment, ie., that 
of the recently derived growth data of the 2dFGRS, 
SDSS-LRG, VVDS and WiggleZ galaxy surveys. Note 
that for the background we use two reference expan- 
sion models namely flat ACDM and DGP respectively. 
The interesting aspect of the latter scenarios is that the 
corresponding functional forms of the Hubble parame- 
ters are affected only by one free parameter, that of the 
dimensionless matter density at the present time Q m o- 
The structure of the article is as follows. In section 2, 
we briefly discuss the background cosmological equa- 
tions. The theoretical elements of the growth index are 
presented in section 3 in which we extend the original 
Polarski & Gannouji (2008) method for a large family 
of 7(2) parametrizations. In section 4 we briefly discuss 
the growth data. In section 5, we perform a likelihood 
analysis in order to constrain the growth index model 
free parameters. Finally, the main conclusions are sum- 
marized in section 6. 



2 THE BACKGROUND EVOLUTION 

For homogeneous and isotropic flat cosmologies, driven 
by non relativistic matter and an exotic fluid (DE mod- 
els) with equation of state (hereafter EoS), pde = 
w(a)poE, the first Friedmann equation can be written 

as: 



H 2 {a) 
HI 



E\a) 



fimOfl + fioB.oe 



f* dlny[l + w(y)] ,^ 



where E(a) is the normalized Hubble flow, a(z) — 1/(1 + 
2) is the scale factor of the universe, w(a) is the EoS 
parameter, Q. m Q is the dimensionless matter density at 
the present time and 0,de,o = 1 — fi m o denotes the DE 
density parameter. Using the Friedmann equations, it is 
straightforward to write the EoS parameter in terms of 
E{a) = H(a)/H (Saini et al. 2000; Huterer & Turner 
2001) 



w(a) = 



where 



O m (a) 



2 dlnE 

3 da 



1 - Q. m (a) 



n, 



:(>(! 



E 2 {a) 



(2) 



(3) 



Differentiating the latter and taking into account eq. (2) 
we obtain 

dQri 



- -w(a)fl m (a) [1 - Qm(a)] 
da a 



(4) 



Since the exact nature of the DE is unknown, the above 
DE EoS parameter includes our ignorance regarding the 
physical mechanism powering the late time cosmic ac- 
celeration. It is also worth noticing that the concor- 
dance ACDM cosmology is described by a DE model 
with w(a) = —1. 

Interestingly, the above method can be generalized 
to the context of modified gravity. Indeed, instead of 
using the exact Hubble flow through a modification of 
the Friedmann equation one may consider an equiva- 
lent Hubble flow somewhat mimicking cq. (1). The in- 
gredient here is that the accelerating expansion can be 
attributed to a kind of "geometrical" DE contribution. 
Now, due to the fact that the matter density (bary- 
onic+dark) cannot accelerate the cosmic expansion, it 
is fair to utilize the following parametrization (Linder 
& Jenkins 2003; Linder 2004): 

E 2 (a) = ^-^ = n m0 a- 3 + AH 2 . (5) 

It becomes clear that any modification to the Friedmann 
equation of general relativity is included in the last term 
of the above expression. Now using eqs. (2) and (5) one 
can derive the effective ( "geometrical" ) dark energy EoS 
parameter 

1 dlnAH 2 
to(o) = -l-- _ ■ (6) 



dlna 

In the context of a flat DGP cosmological model 
the "accelerated" expansion of the universe can be ex- 
plained by a modification of the gravitational interac- 
tion in which gravity itself becomes weak at very large 
distances (close to the Hubble scale) due to the fact 
that our four dimensional brane survives into an extra 
dimensional manifold (see Deffayet, Dvali & Cabadadze 
2002 and references therein). An interesting feature of 
this pattern is that the corresponding functional form 
of the normalized Hubble function as given by eq. (5), 
contains only one free parameter, il m o- The quantity 
AH 2 is given by 



AH 2 



2Q b 



2m()ffl 3 + &b 



(7) 



where Qt, w — (1 — Q m o) /4. From eq.(6), it is readily 
checked that the geometrical (effective) DE equation of 
state parameter reduces to 

1 



w(a) = 



1 + f2 m (a) 



(8) 



In this model due to its gravity nature, the effective 
Newton's parameter G e g is not any more the usual con- 
stant Gn but it takes the following form (Lue, Scossi- 
marro & Starkman 2004) 
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Geff(o) = GjvQ(a) 



Q(a) 



2 + 4gjUo) 

3 + 3Q m (a) 



(9) 



3 THE LINEAR GROWTH RATE 

For the purpose of the present study, we first discuss the 
basic equation which governs the evolution of the matter 
perturbations within the framework of any DE model 
(scalar or geometrical). An important ingredient in this 
analysis is the fact that at the sub-Hubble scales the 
DE component is expected to be smooth and thus one 
can use perturbations only on the matter component of 
the cosmic fluid (Dave, Caldwell & Steinhardt 2002) . In 
particular, following the notations of Lue et al. (2004), 
Linder (2005), Stabenau & Jain (2006), Uzan (2007), 
Linder & Cahn (2007), Tsujikawa, Uddin & Tavakol 
(2008) and Dent, Dutta & Perivolaropoulos (2009) we 
can derive the well known scale independent equation 
of the linear matter overdensity 5 m = Sp m /pm 



6m + 2HS„ 



47rGcffPm<5„ 



(10) 



a solution of which is 5 m (t) oc D(i), with D(t) denot- 
ing the linear growing mode (usually scaled to unity at 
the present time). Notice, that p m is the matter density. 
Of course, for the scalar field DE models [G c a = Gn, 
Q(a) = 1], the above equation reduces to the usual time 
evolution equation for the mass density contrast (Pee- 
bles 1993), while in the case of modified gravity models 
(see Lue et al. 2004; Linder 2004; Linder & Cahn 2007; 
Tsujikawa et al. 2008; Gannouji et al. 2009) we have 
G c fi t^ Gn (or Q(a) 7^ 1). Transforming equation (10) 
from t to a (4^ = H d f n - ) , we simply derive the evolu- 
tion equation of the growth factor D(a) 



a 1 d 2 D / d\nE\ a dD 

+ 3 + a- 



D da 2 



da J D da 



-fi m (a)Q(a) . (11) 



We would like to remind the reader here that solving 
eq.(ll) for the concordance A cosmology , we derive 
the well known perturbation growth factor (see Peebles 
1993) 



D(z) 



5fl m oE(z) 



(1 + u)du 
E 3 (u) 



(12) 



We would like to stress that for the ACDM cosmological 
model we use the above equation normalized to unity 
at the present time. 



3.1 The evolution of the growth index 

As we have mentioned in the introduction, for any type 
of DE, an efficient parametrization of the matter pertur- 
bations is based on the growth rate of clustering origi- 
nally introduced by Peebles (1993). This is 

dlnD 



/W din a 
which implies 



n m (a) 



(13) 



For the usual ACDM cosmological model we have w(a) 
-1, Q^i 11 ) = 1 ~~ ^m(a) and Q(a) = 1. 



D(a) — exp 



i£ m (,-£/ 



dx 



(14) 



where 7 is the so called growth index (see Silveira 
& Waga 1994; Wang & Steinhardt 1998; Lue et al. 
2004; Linder 2004; Linder & Cahn 2007; Nesseris & 
Perivolaropoulos 2008). 

Combining eq.(13), eq.(ll) and eq.(2), we find after 
some simple algebra 



«£ + f + X(a)f = ^n m (a)Q(a) 
where 



X(a) 



-w(a) [1 - n m (o)] 



(15) 



(16) 



If we change variables in eq.(15) from a to redshift [-4- — 
— (1 + z)~ 2 -$-] and utilizing eqs.(4) (13), then we can 
derive the evolution equation of the growth index 7 = 
7(2) [see also Polarski & Gannouji 2008] 



-(1 + 2) 7 'ln(n m ) + W m + 3w(l - n m )(7 -\) + \ 



(17) 



Evaluate eq.(17) at z = we have 

-7'(o)hi(Ono) + affl + 3w (i - n m0 )[7(o) -\] + \ 

= |Qo^7 (0) , (is) 

where Qo = Q(z = 0) and wo = w(z = 0). 

It is interesting to mention here that the last few 
years there have been many theoretical speculations re- 
garding the functional form of the growth index and 
indeed various candidates have been proposed in the 
literature. In this work, we decide to phenomenologi- 
cally treat the functional form of the growth index 7(2) 
as follows 



7(2) =70+712/(2) 



(19) 



In other words, the above equation can be viewed as a 
first order Taylor expansion around some cosmological 
quantity such as a(z), z and Q m {z). Interestingly, for 
those y(z) functions which satisfy y(0) — [or 7(0) = 
70] one can write the parameter 71 in terms of 70. In 
this case [7'(0) = 7iJ/'(0)], using eq.(18) we obtain 



7i 



^°o+3wo(7o-|)(l-«mo) 



3 n o 1_ to _L 1 



y'(0)lnQ mO 
In brief, we present various forms of 7(2), Vz. 



.(20) 



• Constant growth index (hereafter To model): Here 
we set 71 strictly equal to zero which implies 7 = 70. 

• Expansion around z = (see Polarski et al. 2008; 
hereafter Fi model): In this case we have y(z) = z. 
However, the latter parametrization is valid at relatively 
low redshifts < z < 0.5. In the statistical analysis 
presented below we use a constant growth index namely 
7 = 7o + O.571 for z > 0.5. 

• Interpolated parametrization (hereafter F2 model): 
Owing to the fact that the Fi model is valid at low 
redshifts we propose to use a new formula y{z) = zcT z 
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Table 1. The growth data. The correspondence of the columns is as follows: index, redshift, observed growth rate and references. 
In the final column one can find various symbols of the data appearing in Fig. 2. 



Index 



Rcfs. 



Symbols 



1 


0.17 


0.510 ±0.060 


Song 


2 


0.35 


0.440 ± 0.050 


Song i 


3 


0.77 


0.490 ±0.180 


Song 


1 


0.25 


0.351 ±0.058 




5 


0.37 


0.460 ± 0.038 




6 


0.22 


0.420 ± 0.070 




7 


0.41 


0.450 ± 0.040 




8 


0.60 


0.430 ± 0.040 




9 


0.78 


0.380 ± 0.040 





Song & Percival 2009; Percival et al. 2004 open circles 

Percival 2009; Tegmark et al. 2006 open circles 

k Percival 2009; Guzzo et al. 2008 open circles 

Samushia et al. 2012 open triangles 

Samushia et al. 2012 open triangles 

Blake et al. 2011 solid circles 

Blake et al. 2011 solid circles 

Blake et al. 2011 solid circles 

Blake et al. 2011 solid circles 



which connects smoothly low and high-redshifts ranges. 
The latter y{z) formula can be seen as a combination of 
Ti model with that of Dossett et al. (2010). Obviously, 
at large redshifts z>lwe have 700 ~ 70. 

• Expansion around a — 1 (Wu, Yu & Fu 2009; Ishak 
& Dosset 2009; Belloso et al. 2011; Di Porto et al. 2012 
hereafter T3 model): Here we use y(z) — 1 — a(z) — -^j— 
which implies that for z>lwe get 700 — 70 + 71 ■ 

• Expansion around Q. m — 1 (Wang & Steinhardt 
1998 hereafter T4 model): Now we parametrize y(z) as 
follows y(z) — 1 — Q m (z). For the DE models with a con- 
stant EoS parameter w(z) = too one can write (70,71) 
only in terms of wo 

3(1 -too) _ 3 (1 - w )(l - 3wo/2) 



7o 



6to 



7i 



125 



(1 - 6w /5) 3 



21 



At large redshifts fl m (z) ~ 1 we get 700 ~ 70. Note, 
that the DPG cosmological model predicts (70,71) — 



( 



(Linder 2004; Linder et al. 2007; Gong 2008). 



From the above presentation it becomes evident 
that for the Ti_3 parametrizations we have j/(0) = 
and j/'(0) = 1, respectively. Therefore, for the case of the 
ACDM cosmology with (fi m0 ,7o) = (0.273, £) cq.(20) 
provides 71 ~ —0.0478, while for the case of the T4 
model we obtain 71 ~ 0.01127 (see eq.21). In addition, 
based on the DGP gravity with (fi m o,7o) = (0.273, y|) 
the Fi_3 models give 71 ~ 0.05. 



4 THE GROWTH DATA 

The growth data that we will use in this work based 
on 2dF, SDSS and WiggleZ galaxy surveys, for which 
their combination parameter of the growth rate of struc- 
ture, f(z), and the redshift-dependent rms fluctuations 
of the linear density field, <Js(z), is available as a func- 
tion of redshift, f(z)ag(z). The fag = A estimator is 
almost a model-independent way of expressing the ob- 
served growth history of the universe (Song & Percival 
2009). In particular, we will use: 

• The 2dF (Percival et al. 2004), SDSS-LRG 
(Tegmark et al. 2006) and VVDS (Guzzo et al. 2008) 
based growth results as collected by Song & Percival 
(2009). This sample contains 3 entries. 

• The SDSS (DR7) results (2 entries) of Samushia et 
al. (2012) based on spectroscopic data of ~106000 LRGs 
in the redshift bin 0.16 < z < 0.44. 




Figure 1. Left Panel: The variance A\ = X xL,in arounc l 
the best fit 7 value for the A cosmology. Note that the cross 
corresponds to (7Ai^Xi CT ) = (tt>1)- Rvht Panel: The sta- 
tistical results in the case of the DGP model. The corre- 
sponding cross is ("I dgp, AxU = (t1>9). 



• The WiggleZ results of Blake et al. (2011) based on 
spectroscopic data of ~ 152000 galaxies in the redshift 
bin 0.1 < 2 < 0.9. This dataset contains 4 entries. 

In Table 1 we list the precise numerical values of the 
data points with the corresponding errors bars. 



5 FITTING MODELS TO THE DATA 

In order to quantify the free parameters of the growth 
index we perform a standard \ 2 minimization pro- 
cedure between N — 9 growth data measurements, 
A obs = fobs{z)cr&,obs{z), with the growth values pre- 
dicted by the models at the corresponding redshifts, 
Mp> z ) = /(p, z)a% (p,z) with a 8 (p,z) = cr 8 ,oD(j>,z). 
The vector p contains the free parameters of the model 
and depending on the model. In particular, the essential 
free parameters that enter in the theoretical expectation 
of are: p = (70, 71, Umo)- The \ 2 function' is defined as: 



T Likelihoods are normalized to their maximum values. In 
the present analysis we always report lcr uncertainties on the 
fitted parameters. Note also that the total number of data 
points used here is TV = 9, while the associated degrees of 
freedom is: dof= N — k — 1, where k is the model-dependent 
number of fitted parameters. The uncertainty of the fitted 
parameters will be estimated, in the case of more than one 
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N 

X 2 (zi\p) = ^2 

i = l 



A obs (zi) - A(p,z x ) 



(22) 



where at is the observed growth rate uncertainty. To 
this end we will use, the relevant to our case, cor- 
rected Akaike information criterion for small sample size 
(AIC C ; Akaike 1974, Sugiura 1978), defined, for the case 
of Gaussian errors, as: 



AIC C 



2 „, 2fc(fc 

Xmin "T 2/C + 



N-k- 



(23) 



where k is the number of free parameters, and thus when 



k = 1 then AIC C 



Xr> 



+ 2. A smaller value of AIC C 



indicates a better model-data fit. However, small differ- 
ences in AICc are not necessarily significant and there- 
fore, in order to assess, the effectiveness of the different 
models in reproducing the data, one has to investigate 
the model pair difference AAIC C = AIC C ,;, — AIC C , X . 
The higher the value of |AAIC C |, the higher the ev- 
idence against the model with higher value of AIC C , 
with a difference | AAIC C |> 2 indicating a positive such 
evidence and |AAIC C |<; 6 indicating a strong such ev- 
idence, while a value <, 2 indicates consistency among 
the two comparison models. A numerical summary of 
the statistical analysis for various 7(2) parametrizations 
is shown in Table 2. In general, we find that our results 
are in agreement, within la uncertainties, with previ- 
ous studies (Di Porto et al. 2008; Gong 2008; Nesseris 

6 Perivolaropoulos; Dosset et al. 2010; Fu et al. 2009; 
Basilakos 2012). 

5.1 Constant growth index 

First of all we consider the To parametrization (7 = 70 , 
71 = 0; see section 3.1) which implies that the corre- 
sponding statistical vector becomes: p = (7, 0, fl m o). 
We will restrict our present analysis to the choice 
(Qmo, 0-8,0) = (0.273,0.811) provided by WMAP7 (Ko- 
matsu et al. 2011)''. Note that we sample 7 £ [0.1, 1.3] 
in steps of 0.001. 

In the left panel of Fig. 1 we show the varia- 
tion of Ax 2 = X 2 (l) ~ Xminil) around the best 7 fit 
value for the concordance A cosmology. We find that 
the likelihood function of the growth data peaks at 

7 = 0.602 ± 0.055 with xLn ^7.1 for 7 degrees of 
freedom^. Alternatively, considering the ACDM theo- 
retical value of 7 (= 6/11) and minimizing with respect 
to Qmo we find i} m o = 0.243 ± 0.034 (see also Nesseris 

6 Perivolaropoulos 2008) with xLin/dof ~ 7.37/7. 
Our growth index results are in agreement within la 
errors, to those of Samushia et al. (2012) who found 

7 = 0.584 ± 0.112. However, our best-fit value is some- 
what greater and almost la (Axi a — 1) away, from the 
theoretically predicted value of 7A -^ 6/11 (see cross 



such parameters, by marginalizing one with respect to the 

others. 

t For the DGP model Gong (2008) found Q m0 = 0.278. 

8 Using eq. (14) in the likelihood analysis for the usual A 

cosmology we obtain 7 = 0.595 ± 0.071 with x 2 /d°f — 

7.59/7. Note that for the DGP model we only use eq.(14). 



_ACDM 
_DGP 




£*-*^ 



Figure 2. Comparison of the observed and theoretical evo- 
lution of the growth rate A(z) = f(z)ag(z). The solid and 
dashed lines correspond to ACDM (7 = 0.602) and DGP 
(7 = 0.503) expansion models respectively. The thin-line er- 
ror bars correspond to la 7-uncertainties for the A cosmol- 
ogy. We do not plot the la 7-uncertainties for the DGP model 
in order to avoid confusion. The different growth datasets are 
represented by different symbols (see Table 1 for definitions). 



in the left panel of Fig. 1). It is interesting to mention 
here that such a small discrepancy between the theo- 
retical ACDM and observationally fitted value of 7 has 
also been found by other authors. For example, Di Porto 
et al. (2008) obtained 7 = 0.60±g;|g, Gong (2008) mea- 
sured 7 = 0.64^0^5 while Nesseris & Perivolaropoulos 
(2008) found 7 = 0.67+J]; 2 ?. Recently, Basilakos (2012) 
and Hudson & Turnbull (2012) using a similar analysis 
found 7 = 0.613to;og3 and 7 = 0.619 ± 0.054 respec- 
tively. 

Concerning the DGP model (see the right panel of 
Fig. 1) the best fit parameter is 7 = 0.503 ± 0.06 with 
Xrnin/dof ~ 5.32/7. If we fix the value of 7(= 11/16) to 
that predicted by the DGP model we find a rather large 
value of the dimensionless matter density at the present 
time, n m0 = 0.380 ± 0.042 with xLin/dof ~ 5.38/7. 

The value of A1Cc,dgp(— 7.32) is smaller than 
the corresponding ACDM value which indicates that 
the DGP model (7dgp = 0.503) appears now to fit 
slightly better than the usual A cosmology the growth 
data. However, the small |AAIC C | value (ie., ~ 1.8) in- 
dicates that the two comparison models represent the 
growth data at a statistically equivalent level. On the 
other hand form the right panel of Fig.l, it becomes 
clear that the best-fit 7 value is much lower and almost 



3<7 (AxIct — 9) away, from jdgp — 11/16 (see cross 
in the right panel of Fig. 1) implying that the growth 
data disfavor the DGP gravity. We would like to stress 
here that the above observational DGP constraints are 
in excellent agreement with previous studies. Indeed, 
Wei (2008) found 7 = 0.438JT^ ;^. Also Gong (2008) 
and Dosset et al. (2010) obtained 7 = 0.55+q Jg and 
7 — 0.483^0088 respectively. In Fig. 2, we plot the mea- 
sured A i,s(z) with the estimated growth rate function, 
A(z) = f(z)a s (z) [see ACDM - solid line and DGP - 
dashed line]. 

The goal from the above discussion is to give 
the reader the opportunity to appreciate the relative 
strength and precision of the different methods used in 
order to constrain the growth index. It becomes evident 
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Figure 3. The ACDM expansion model: Likelihood contours 
(for Ax 2 = -21n£/£ max equal to 2.30, 6.18 and 11.83, corre- 
sponding to la, 2cr and 3c confidence levels) in the (70,71)- 
The upper left and right panels show the results based on the 
Ti— 2 paramctrizations. In the bottom left and right panels 
we present the contours of the T3_4 parametrizations (for 
more details see section 3.1). We also include the theoretical 
ACDM (70,71) values given in section 3.1. 



that with the combined high-precision /erg growth rate 
data of Song & Percival (2009), Samushia et al. (2012) 
and Blake et al. (2011) we have achieved to place quite 
stringent constraints on 7. 



5.2 The Ti_4 parametrizations 

After we have presented the simplest version of the 
growth index, it seems appropriate to discuss the obser- 
vational constraints on the time varying growth index, 
j(z). Following the considerations exposed in section 
3.1, hereafter we will set p = (70,71,0.273) in equa- 
tion (22). In Fig. 3 (ACDM model) and Fig. 4 (DGP 
model) we present the results of our statistical analysis 
for the Ti (upper left panel), F-2 (upper right panel), 
Ts (bottom left panel) and F4 (bottom right panel) 
parametrizations in the (70,71) plane in which the cor- 
responding contours are plotted for la, 2a and 3a con- 
fidence levels. Notice, that we sample 70 £ [0.1, 1.3] 
and 71 G [-2.2, 2.2] in steps of 0.001. The theoreti- 
cal (70,71) values (see section 3.1) in the ACDM and 
DGP expansion models indicated by the crosses. Over- 
all, we find that the predicted ACDM (70,71) solutions 
of the Fi_4 parametrizations remain close to the la 
borders (Axi CT — 2.30; see crosses in Fig. 3). Regard- 
ing the DGP model (see Fig. 4) we would like to stress 
that the predicted (70,71) values approach the 3a bor- 
ders (Ax| CT — 11.83; see crosses on Fig.4) of the 70 — 71 
contours. Obviously, this is a clear indication that the 
current growth data can not accommodate the DGP 
gravity model. 

Below we briefly discuss the main statistical re- 
sults: (a) Pi parametrization: For the usual A cosmol- 

0.40lH 8 80 ^d 



ogy the likelihood function peaks at 70 



7i 



0.603 ± 0.241 with Xmin/dof ~ 5.74/6, while in 



Figure 4. The Likelihood contours for the DGP expansion 
model (for more definitions see caption of figure 3). Here 
the crosses correspond to the theoretical DGP (70,71) pair 
provided in section 3.1. 



and 71 = 0.164 ± 0.221 with xtnin/dof ~ 5.10/6. Inter- 
estingly, the use of the combined growth data provides 
a significant improvement in the derived (70 , 71 ) con- 
straints with respect to the previous studies (Di Porto et 
al. 2008; Gong 2008; Nesseris & Perivolaropoulos 2008; 
Dosset et al. 2010; Fu et al. 2009). 

Such an improvement is to be expected because the 
WiggleZ and the SDSS-DR7 surveys measure f(z)a$(z) 
to within 8 — 17% (Blake et al. 2011; Samushia et al. 
2012) in every redshift bin, in contrast to the old growth 
rate data (Song & Percival 2009) in which the corre- 
sponding accuracy lies in the interval 12 — 37%. 
(b) Now we concentrate on the F-2 and T3: parametriza- 
tions: We find that within la errors we can put some 
constraints on the free parameters. In particular, the 
best fit values are: (i) ACDM: for T2 we have 70 = 
0.311±g:g|g, 71 = 1.221 ±0.343 (xlijdof ~ 5.26/6) 
while for F 3 we get 70 = 0.345±g;g|g , 71 = 1.006 ±0.314 
{Xmin/dof ~ 5.06/6) and (ii) DGP: in the case of 



T2 model we obtain 70 = 0.401 



-0.094 r 
-0.090' 71 



0.384 ± 



0.320 (xLin/dof ~ 5.00/6) and for T 3 we find 70 = 
0.412+° ° f , 71 = 0.321 ± 0.290. In the latter case the 
reduced xL» is ^ 4.94/6. 

(c) T4 parametrization: In this case the 70 is strongly 
degenerate with 71 (see bottom left panels in Figs. 3, 4). 
Indeed, we can provide the following 71 — 70 relations: 



7i = 



1.134(±0.005) - 1.879(±0.006) 7 o ACDM 
0.887(±0.005) - 1.749(±0.006)7o DGP 



(24) 



the case of the DGP gravity we obtain 70 = 0.441_ 



Finally, as we have already mentioned in Table 2, 
one may see a more compact presentation of our statis- 
tical results. For both cosmological (ACDM and DGP) 
models the information theory pair model characteri- 
zation parameter, AAIC C , indicates that all the 7(2:) 
functional forms explored in this study are statistically 
equivalent in representing the growth rate data, since 
JAAIC C | < 2 for any pametrization pair. In Fig. 5 we 
present the evolution of the growth index for various 
parametrizations. In the case of the concordance A cos- 
mology (upper panel of fig. 5) the relative growth in- 
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Table 2. Statistical results for the combined growth data (see Table 1): The I s * column indicates the expansion model, the 2 nd 
column corresponds to 7(2) paramctrizations appearing in section 3.1. 3 rd and 4 rth columns show the 70 and 71 best values. 
The remaining columns present the goodness-of-fit statistics (reduced \ 2 an d AIC C ). 
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Figure 5. Upper Panel: The evolution of the growth index 
for the ACDM model. The lines correspond to To (solid), 
Ti (short dashed), r 2 (long dashed) and T3 (dot-dashed). 
The error bars correspond to lcr 7-uncertainties for the r 2 
paramctrization. We do not plot the la 7-unccrtainties for 
the other paramctrizations in order to avoid confusion. Bot- 
tom Panel: The evolution of the growth index for the DGP 
cosmological model. 



dex difference of the various fitted 7(2) models indicates 
that the Fi _3 models have a very similar redshift depen- 
dence for z < 0.5, while the F3 parametrization shows 
very large such deviations for z > 0.5. Based on the 
DGP gravity model (bottom panel of fig.5) we observe 
that the ri_ 2 parametrizations provide a similar evolu- 
tion of the growth index. The T3 paramctrization shows 
large deviations at large redshifts z > 1.5. However the 
large 7(2) errors appear in fig.5 are due to the large un- 
certainty of the 71 fitted parameter, implying that more 
and accurate data are essential in order to distinguish 
among the different 7(2) functional forms. 



6 CONCLUSIONS 

It is well known that the so called growth index 7 plays 
a key role in cosmological studies because it can be 
used as a useful tool in order to test Einstein's gen- 
eral relativity on cosmological scales. We have utilized 



the recent growth rate data provided by the 2dFGRS, 
SDSS-LRG, VVDS and WiggleZ galaxy surveys, in or- 
der to constrain the growth index. Performing a likeli- 
hood analysis for various 7(2) parametrizations, we ar- 
gue that the use of the above combined growth data 
places the most stringent constraints on the value of 
the growth index. Overall, considering a ACDM ex- 
pansion model we find that the observed growth in- 
dex is in agreement, within lcr errors, with the theo- 
retically predicted value of 7a ~ 6/11. In contrast, for 
the DGP expansion model we find that the measured 
growth index is almost 3<r away from the corresponding 
theoretical value ^dgp — 11/16 which implies that the 
present growth data can not accomodate the DGP grav- 
ity model. Finally, considering a time varying growth in- 
dex parametrization namely 7(2) = 70 + 712/(2) [where 
y(z) — z,ze~ z ,l — a(z) and 1 — Q, m {z)] we find that 
although the 70 parameter is tightly constrained, the 71 
parameter remains weakly constrained. Hopefully, with 
the next generation of surveys, based on Euclid and Big- 
BOSS, we will be able to put strong constraints on 71 
and thus to check departures from 7 = const. 
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